a \ Middle East Journal of Applied Science & Technology (MEJAST) 
M E J A S T Volume 6, Issue 1, Pages 38-41, January-March 2023 


Shilpa kulkarni! & Pralahad Mahagaonkar” 


‘Assistant Professor, Department of Physics, SSA Govt First Grade College (Autonomous), Ballari, Karnataka, India. = 
?Assistant Professor, Department of Mathematics, BITM, Ballari, Karnataka, India. es Crossref 
Corresponding Author (Pralahad Mahagaonkar) - pralaahdm74@ gmail.com* 


DOI: https://doi.org/10.4643 1/MEJAST.2023.6104 


Copyright © 2023 Shilpa Kulkarni & Pralahad Mahagaonkar. This is an open-access article distributed under the terms of the Creative Commons Attribution 
License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 


ABSTRACT 


. Introduction 


Now a days it has been observed in the science and technology that the development of computational methods are 
a part of research area to solve the linear differential equations. The analysis was made in the applied sciences and 
few of them solved related to numerical methods in the field of science. The use of numerical methods have seen to 
draw a great attention in science and engineering field. Many of the researchers published their research papers in 
these fields. Many methods are there to solve damped and undamped force of vibrations. In this paper undamped 
force vibration of a spring have been solved by linear differential equations by using analytical and numerical 


methods. 


Undamped force of vibration defined as “When there is no external force act on a vibrating body, if it is 
continuously vibrating without damping then this system of vibration is called undamped force of vibration’’. In 
undamped force vibration there is no energy dissipation to the surrounding system. It is a simple harmonic motion 
(SHM) oscillation with constant amplitude with respect to time. Examples alternating current and alternate voltage, 


swinging pendulum in vacuum. 


“2. Undamped force vibration of spring 


Consider the undamped force vibrations of spring given by the differential equation is, 


d°x 
Mg PG) (1) 


In this paper we take the special choice of f(t)=(1-sint), m=1kg, k=1N/m, with initial conditions x(O) = x'(0) =0 
then equation (1) gives us, 


2 
T+ x(t)=(1-sin2), 


As the initial conditions, 
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The exact solution of equation (1) by using the classical method is, 
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ti .. md # 
x(t) =—+sint -——+— 
2 2 2 


Applying initial conditions we get C, =—1 , C, =O. 

3. Runge-Kutta Fourth Order Method 

To compute y(x, + A) for the required value of y'(x) +h) = z(x) +h). 
Firstly, we compute the values of 

Kk, = hf (%; Yo» 2) : 1, =hg (Xo, Yo» 2) 
ky =hf (x) +h12, yo +k, 12,2 +1, /2) 5 Lb =hg(x) +A/2, yy) +k, 12, 2% +1, /2) 
ks = hf (x) +h/2, yo + ky 12, 2% +1,/2) > b= helm th/2,y,4+,/2,%, 41/2) 
ky =f (% +s Yo + hq2% +4) pba He hy Ee Va thy Zy FL) 
y(X% +h) = yy +1/6(k, + 2k, + 2k, +k,) and 

y (X) +h) = 2(x, +h) = z, +1/6(L, + 2, + 21, +1,) 


“4. Eulers Modified Method 


If h is the step size, Let’s consider equation differential equation of (1) 


x =(1—-sint)—x(ft) (2) 

dx : 

= =z=f(t,z,x) then equation (1) reduces to, 
& _ a-sin())—x= g(t,z,x) 

dt 


Consider the Eulers approximation values as follows, 
X, =X +hf (to, Zos%) 
Zy = Zqy FAG (bys ZsXp) aNd SO OM... 


Table 1. Results from Different numerical methods 


h=1/32 Exact Solution | Runge—Kutta Method x Eulers Method 
t T Error t Error 
1 1.2835 1.30891 0.42758 1.3005 0.43128 
3 1.2568 1.74521 0.52969 1.72412 0.52712 
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a) 2.0012 2.18151 0.53833 2.1717 0.52181 
7 3.1258 3.05411 0.43125 3.0521 0.42181 
9 2.8592 3.49041 0.47056 3.50125 0.42158 
11 3.2561 4.36301 1.06535 4.32156 1.06421 
13 3.1258 4.79932 1.6196 4.80012 1.62131 
5 4.1258 5.23562 2.24149 5.28523 2.25863 


5. Results and Discussions 


In the above table 1 results have been analyzed exact solution method with the numerical methods. For both the 
numerical methods the step size was fixed as h = 1/32. It has been noticed in the results that whenever step size was 
increased the estimated error has also been increased. The studied results showed in the above listed table 1 such 


that the time period value was different for different methods. 
“6. Conclusion 


Here we have assumed and taken the equation of undamped force of vibrations in the form of differential equation 
and hence solved analytical method and numerical method. The comparison between these methods have been 


discussed. 
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